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We consider the small-angle multiple neutron scattering and a possibility of its model-free anal- 
ysis by the inverse problem method. We show that the ill-defined inverse problem is essentially 
regularized by the use of a planar detector of neutrons without a "beam-stop" facility. We analyze 
the obtained expressions with the examples of different models of scatterers. 
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The small angle neutron scattering (SANS) technique 
is an important tool for studying large-scale inhomo- 
geneities in condensed matter. The profile of the in- 
tensity of passing neutrons around the (0, 0, 0) reflection 
provides an information about the size and shape of the 
inhomogeneities. As discussed, e.g., in [[[], the procedure 
of extracting this useful information may be complicated 
by different factors, the depletion of the incident beam 
and multiple scattering among them. 

Different ways of analysis of the multiple scattering 
were proposed up to now j2|^,|l[, however all of them 
depend to some extent on the particular assumptions. 
In fact, the assumptions are almost inevitable, since the 
analysis of multiple SANS data deals essentially with an 
inverse problem, which is mathematically ill-defined. 

In this paper we show that the degree of uncertainty 
in the analysis of multiple scattering can be substantially 
reduced, when one uses the neutron detector without a 
so-called "beam-stop" facility. This latter item is imple- 
mented in many experimental setups in order to decrease 
the presumably useless contribution of neutrons passed 
without a collision. We argue that when a notable part 
of neutrons scatter in the specimen a few times, this usu- 
ally discarded intensity can provide a crucial ingredient 
for the model-free data analysis. Mathematically, taking 
into account the intensity at the zero angle corresponds 
to the substantial regularization of the inverse problem 
found in the multiple SANS. 

The rest of this paper is organized as follows. First we 
rederive the basic equation known in the Bethe-Moliere 
theory of multiple scattering. Our treatment here is simi- 
lar to one in jlj . We show an intermediate formula which 
is more suitable for the desired inversion. Formally in- 
verting this formula, we discuss the obtained expression 
to some detail. The form of functions expected in some 
models of scattering is shown and discussed in the last 
part of the paper. 

The scattering probability is determined by the sin- 
gle scattering cross section a = da/dfl with thee mo- 



mentum transfer Q. Given the momentum of the in- 
cident neutron hi, it is convenient to define the trans- 
verse "momentum" q = Q/fc; and the normalized quan- 
tity cr(q) = cr(q)/(J d 2 qa(q)). The value c(q) is a 
conditional probability of the scattering at the angle q. 
The total scattering cross section erg = j d 2 qer(q) en- 
ters the definition of a mean- free path, /, in a usual way, 
/ = (aorii)^ 1 , where rii is the scatterers' density. 

In the diffractional limit, when I is much larger than 
the typical size of inhomogeneities, the probability of 
multiple scattering is given by a Poisson distribution 
so that n scattering events on the path L are ex- 
pected for the fraction p n of incoming neutrons, p n = 
(L/l) n exp(— L/l). The partpo = exp(—L/l) corresponds 
to the non-scattered neutrons. 

The total multiple scattering intensity is given by a 
sum 



oo 
n=0 



with the recursively defined I n (q) : 
Jo(q)=*(q) 

4(q) = / d 2 k/, i _ 1 (q-k) ( r(k), 

In terms of 2D Fourier transform A(r) 
we rewrite 



J„(r) = [a(r)] n , n >0 



The sum (|J) is explicitly represented in the form 
I(r,L) = exp(-L/l)J2(jv(r)) 



exp 



(l-a(r)) 



(1) 
(2) 

n > 1 (3) 
= J dqe ic *M(q) 

(4) 
(5) 

(6) 
(7) 
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Note that, according to our definition, er(r = 0) = 1. 
Let us assume now that the scattering is isotropic, i.e., 
cr(q) does not depend on the direction of q. It follows 
then that the quantity tr(r) depends only on the absolute 
value of r, and we can transform (Q) by going to the "mo- 
mentum" representation and integrating over the angle 
in r— plane. As a result, we obtain 



f°° rdr 
i (q, l ) = / TT- Jo(qr)exp 
Jo 27r 

= S(q)e- L / 1 - 



-y(l-a(r)) 



* dr 
2V 



Jo(qr) 



e *(r)L/l _ l 



,-L/l 



with the Bessel function Jo(x). In the second line of (j8|) 
the integral is regularized by explicitly subtracting the 
singular contribution of the direct beam. One easily ver- 
ify the equivalence of the expression (@) to the Bethe's 
formula for the small-angle scattering. [Q 

The Eq. (||) has two disadvantages. First, it demands 
the isotropy of the scattering in the xy— plane, which 
could be violated, e.g., for the textured compounds. Sec- 
ond, Eq. (|^) is rather awkward, which makes hard its 
subsequent analysis with taking into the account the ex- 
perimental resolution function. This latter shortcoming 
of (|^) cannot be cured if one works experimentally with 
the stripe (one-dimensional) neutron dectector which was 
a common situation a few decades ago. On the other 
hand, the most of the neutron centers are equipped nowa- 
days by two-dimensional (2D) SANS machines, and it 
provides a possibility for the easier data analysis as ex- 
plained below. 

Keeping in mind the 2D experimental geometry, we use 
the Eq. ([?]), rather than its particular form (||). Next 
we consider the experimental situation when the ideal- 
ized intensity /(q, L) is convoluted with the instrumen- 
tal resolution function -F(k). For simplicity, we assume 
that this is the Fourier convolution. In addition, we take 
into account that some part of neutrons could be lost 
due to absorption or incoherent scattering. These effects 
are irrelevant to the discussed small-angle scattering and 
may be modelled by the unique absorption length, l a . 
The measured (experimental) intensity I exp (q, L) is then 
given by 



p (q, L) = e- L ^ J d 2 kl(q - k, L)F(k), 



-L/U 



I(r,L)F(r). 



(10) 



(11) 



At L — one has I(r) = 1 and the obvious equality 
I exp (r,L — 0) = F(r), i.e. without a sample one mea- 
sures only the shape of instrumental resolution. 

It should be stressed here, that the next step of our 
analysis will require a certain conservation rule. Specif- 
ically, we observe that in the absence of absorpition 
(L/l a -C 0), one has I exp (r = 0, L) = 1, which equality 



means simply that the total number of neutrons, scat- 
tered and passed without collision, does not depend on 
the thickness L. This obvious conservation law is vi- 
olated, when one excludes from the consideration the 
forward-going neutrons. In practice, it is done by the con- 
ventional setup with a "beam-stop" , which assumes par- 
ticularly that the forward-going neutrons carry no useful 
information. We argue here that the part I exp (q ~ 0, L) 

/ s is very important, enabling one to analyze the multiple 
SANS without further adjustable parameters. 

Indeed, given our assumptions above, we see that if the 

(9quantity I exp (r,L) is determined for two samples with 
thicknesses L%, L2, then the ratio I exp (r, Li)/ I exp (r, L2) 
does not depend on F(r). It allows us to invert the equa- 
tions (0), ( |ll| ) and to write 



1 



Iexp{?, Li) 

L2 — Li I exp (r, L2) 



log 



1 - <r(r) 
I 



/(r) (12) 



This equation is the main result of our paper. Let us 
discuss it in more detail. 

1) The rhs of Eq.([l2]) does not depend on the sample 
thickness Li,2, therefore one can combine the scattering 
data, obtained for different Li and to improve the deter- 
mination of the shape of /(r). 

2) For the case of the isotropic scattering, cr(r) = <r(r), 
one averages over the directions in xy— plane in ( |l2| ) and 
not in (||). 

3) Generally, I exp (r,L) is a complex- valued function. 
However, the ratio I exp (r,L)/I exp (r,0) entering ( |l2] ) 
should be real, provided cr(k) = cr(— k). This con- 
dition means simply the inversion symmetry which is 
generally present on the SANS scale (violating proba- 
bly in some special cases, e.g. for some substances in 
the magnetic field). Adopting that er(k) = <r(— k), we 
note that the imaginary part of logarithm in ( |l2[ ) can 
stem from two reasons. The first reason is the sta- 
tistical noise and the second one is the displacement 
of the direct beam (q = 0) during the measurements 
with different Lj. If qdisp is the displacement of the di- 
rect beam, then I exp (r, L)/I exp (r, 0) acquires the factor 
exp(zqdi sp r). This factor can be easily determined by the 
appropriate analysis of the data. The remaining com- 
plexity in the value of the discussed ratio is due to the 
statistical noise. Both discussed sources of the complex 
quantity in the lhs of (|l^) are irrelevant to the physical 
quantity cr(r) and can be stripped off by writing 

fregi?) = T 1 r log 



L2-L1 



lexpi?, L2) 



(13) 



This definition of f reg (r) is the partial regularization of 
our inverse problem, since it projects the experimental 
data into a definite (symmetrical) class of the functions 
cr(q). Hereafter we assume that /(r) is regularized ac- 
cording to (|l3|) and drop the subscript in the function 
f reg (r). In fact, the regularization ( |l3[ ) is not very re- 
strictive. The crucial regularization of the problem ap- 
peared in (|Tl]), corresponding to the assumption of the 
planar detector without " beam-stop" . 
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4) Two basic features in the determined shape of f(r) 
are of interest. First, since cr(0) = 1 one has /(0) = l/l a , 
i.e. the absorption length is found in the same experi- 
ment as the other quantities. Second, one expects that 
cj(oo) = (see below) so that /(oo) = l/l + l/l a and 
hence the value of I is determined. 

5) Knowing I and l a , we find er(r) and can compare the 
form of this function with the different models of scatter- 
ers, partly discussed below. The characteristic scale ro 
of the variation of <r(r) is connected with the size of the 
inhomogcncitics R by a simple relation R — r /ki. 

6) The above limit ct(oo) = corresponds to the van- 
ishing angle between the neighboring neutron detectors, 
which is evidently never realized. In practice, one deals 
with the finite Fourier transform, and the finite-size ef- 
fects appear at largest r, determined by the resolution 
Aq. One may judge confidently about the limiting value 
/(oo) when the curve f(r) saturates at r < l/Aq. This 
saturation takes place when the average size of the scat- 
terers is smaller than l/Aq. In this case one uniquely 
determines two main parameters / and r , along with the 
whole r— dependence of cr(r). 

7) It was proposed in [H] to analyze the multiple scat- 
tering by inverting Eq. ^ with the subtracted S— function. 
We note here that such procedure requires the knowledge 
of the ratio L/l and should be sensitive to its choice. The 
inversion of (0) is also vulnerable due to the ignorance of 
the absorption, l a , and the experimental resolution. 

It is worth to provide here some examples of the possi- 
ble shapes of cr(r). We do it first for the scattering on the 
spheres and for the phenomenological "Ornstein-Zernike- 
squared" cross-section. 

In the first case, one has the spherical scatterers of one 
radius ro , and the differential cross-section in the diffrac- 
tion limit reads as 

osphio) « <r 3 ^3/2(<? r o) (14) 

For this function one finds the image <r(r) in the form 

<r sp h(r) = g sp h{r 2 /irl) (15) 

with H 

g sp h(x) = Vi-x (l + |j 

-» (>-!)* 



i + vT^af 



(16) 



at x < 1 and g sp h{x) = otherwise. 
In the second case, one writes 

<JOS2{q) cx (q 2 rl +9/2) 



(17) 



with the choice of the factor 9/2 explained below. The 
image a(r) is found as 



C0S2M = gos2(r/2r ) 
goS2(x) = {Zxly/2)K x $x/y/2), 



(18) 
(19) 



with the modified Bessel function Ki(x). We note that 
the expression ([l7]) corresponds to the Guinier radius 
R 2 = 4rp/3, which is approximately two times larger 
than in the case of spheres, R 2 g = 3^/5. 

The behavior of the corresponding functions (1 — cr(r)), 
entering eq. (|l^), is shown in Fig. |J. An obvious satura- 
tion at large r takes place in the case of spherical scatter- 
ers, while this feature is not so prominent for the function 
(1 — <Jos2 if))- O n f ne other hand, the point of intersec- 
tion of both curves at (l — a{r)) ~ 0.88 at r ~ 1.4ro turns 
out to be semi- invariant. 

To explain the last statement, we remind that in prac- 
tice one can rarely expect that all the scatterers are of 
the same size. To model the situation, we suggest that 
the sizes of the spherical scatterers are distributed within 
some range, with the probability P(R) to find a sphere 
of radius R. The function P{R) is normalized and the 
mean square of R is r 2 : 



/■OO I* oo 

/ P{R)dR =1, / R 2 P{R)dR = rl 
Jo Jo 



(20) 



Then one can write 



poo 

a{r) = / dRg sph (r 2 /4R 2 )P(R) (21) 
Jo 

The simplest one-parameter family of functions P(R), 
obeying (f2(i|), is given by 



9 / E)2 

P a (R) = -F c ' " 



T[a] 



(22) 



The parameter a characterizes the width of distribution, 
according to a relation (i? 4 ) — (R 2 ) 2 — r^/a, with (. . .) 
standing for the average. Thus the case a — > 00 corre- 
sponds to the previous case of the monodispersed situa- 
tion, Poo{x) = S(R — ro), while P a {x) becomes uniform 
distribution in the limit a — > 0. 

Substituting ( f22| ) into (21), we find the family of curves 
a a {r) in the form 



T[a) 



G 



3/2,3 
0, 1, a 



(23) 



with the Meijer function, implemented, e.g., in Math- 
ematical,. We numerically found that the above point 
(1 — <j(r)) ~ 0.88 and r ~ 1.4ro is semi-invariant for 
the family ( |2^ ) for 0.4 < a < 15. For completeness, we 
also present here the form of cr(q). Up to the overall 
coefficient one has : 



„-2 r 21 
1 32 



q 2 r 2 



0. 3/2,3 

1, a 



(24) 



Let us now discuss our choice of the factor 9/2 in Eq. 
@. This factor appears quite naturally when one poses 
oneself a following question. Which distribution P(s) in 
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( pl| ) leads to the intensity profile of the form (q 2 r 2 + 1)~ 2 
and what is the relation between ro and r here ? 

The answer to this question is based on the observa- 
tion, that cq. ( ^l|) is the Mellin convolution of the func- 
tions g S ph{x) and P(x) and, in principle, may be inverted 
for the known 17(7"). The corresponding derivation in- 
volves the theory of Mellin-Barnes integrals, discussed, 
e.g. in gj. Without going into details, we quote here the 
final result of the form 



P(R) = (2/R)F(R 2 /r 2 ) 

T 



F(x) = 



-3\/25 



1 



(25) 



The behavior of this function is shown in Fig. |2|. There- 
fore the distribution (^5|) of spherical scatterers with av- 
erage (R 2 ) = 7q leads to the intensity profile ([ft]), (|l8|). 

Consider also a situation when the inhomogeneities in 
the specimen form clusters, whose size is described by 
fractal statistics. ||] In the case of simple volume frac- 
tals, the commonly adopted density-density correlation 
function is given by 



(p(r)p(O)) cxr D - 3 exp(-r/0, 



(26) 



with a fractal dimension D and maximum cluster size £. 
Given the function (|26|), one finds the form of the single 
scattering cross section da/dfl, exhibited elsewhere |3|. 
Its counterpart in r— space, entering our eq. (12) reads 
as 



= GS ( w 



1/2 

u i 2 ' 2 



(27) 



The saturation of this function at large r is realized 
when r > £. In particular, it means that one cannot 
uniquely determine the parameters £, I, l a in the case 
when £Ag > 1. 

In conclusion, we discuss the inverse problem for the 
multiple small-angle neutron scattering. We show that 
the consideration of the forward-going neutrons corre- 
sponds to the essential regularization of this problem. 
The proposed formalism may be applied in the region 
where the sample thickness is of order of a few mean-free 
paths I. 
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FIG. 1. Two examples of possible functions in eq. (|l2"|). 
The profiles described by eq. ( |l5[ ) and eq. (|l^) are shown by 
the solid and dashed lines, respectively. 




FIG. 2. The distribution function of the spherical scatter- 
ers, leading to the intensity (|l7]). The function P(r) is mea- 
sured in units of tv" 1 . 
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